The internal Josephson effect in a Fermi gas near a Feshbach resonance 
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We consider a two-component system of Fermi atoms and molecular bosons in the vicinity of a Feshbash 
resonance. We derive an effective action for the system, which has a term describing coherent tunneling of the 
molecular bosons into Cooper pairs and vice versa. In the equilibrium state, global phase coherence may be 
destroyed by thermal or quantum phase fluctuations. In the non-equilibrium regime, the system may show an 
internal AC Josephson effect leading to real time oscillations in the number of molecular bosons. 
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Introduction. Trapped dilute cold Fermi systems 01 are 
one of the most exciting areas of research in modern con- 
densed matter physics. This field off'ers a great variety of 
novel phenomena P,*?, '^j't'j's'i and presents serious chal- 
lenges for both experimentalists and theorists. The magnetic- 
field induced Feshbach resonance 0] provides an unprece- 
dented degree of control over the inter-particle interactions 
as well as the rate at which the interactions are changed. 
By sweeping magnetic field, one can tune the sign and the 
strength of the interaction fS] and experimentally access the 
BCS-BEC crossover physics [10]. Apart from this, by ad- 
justing the sweeping rate, one can drive the system into var- 
ious non-equilibrium states in which different bosonic and 
fermionic species co-exist and interact with each other. 

Interaction between fermions consists of two channels: an 
inelastic channel involving the formation of a molecular state 
of two fermions and a resonant elastic scattering channel. 
Typically, the elastic scattering length is an unremarkable 
function of the applied field. However, in the vicinity of a 
certain value of the field Bq, the scattering length changes 
dramatically 0,|9[l- On the BEC side of the Feshbach reso- 
nance, the molecular energy level is located below the contin- 
uum states and the molecular state is stable. On the BCS side 
of the resonance, high magnetic field breaks up the molecu- 
lar states and one would expect that the molecular component 
vanish. However, recent experimental results |7] and theo- 
retical works suggest that even on the BCS side of 
the resonance, the fraction of molecular states is non-zero. 
One may assume the existence of an energy barrier between 
the molecular state and the continuum and so the molecu- 
lar state on the BCS side of the resonance may actually be 
metastable. This may happen in part due to a stabilizing effect 
of the Fermi liquid coupled to the molecular component. In 
turn, the molecular component affects the fermionic degrees 
of freedom renormalizing effective fermion-fermion interac- 
tions OH- 

At low temperatures, attractive interaction between 
fermions should lead to the formation of Cooper pairs |4, |3> 
The condensate of Cooper pairs would co-exist and in- 
teract with the system of preformed molecular bosons if the 
life-time of the molecular state is long enough. It is inter- 
esting to describe the co-existence of the two coherent states, 
and this issue is the main subject of the present work. In this 



Letter, we study theoretically a Fermi-Bose mixture and dis- 
cuss quasiequilibrium and non-equilibrium properties of the 
system. 

We start with an effective low energy Hamiltonian describ- 
ing the system near the resonance (on the high-field BCS side) 
I151I1(tIi . In this Hamiltonian, the resonant and non-resonant 
processes are separated by introducing an effective Feshbach 
coupling of the fermion component to the metastable molec- 
ular field. We integrate out the one-particle fermionic de- 
grees of freedom and explicitly derive an effective action for 
the system. This action contains a term which can be inter- 
preted as an internal Josepshon tunneling fivll of molecular 
bosons into Cooper pairs and vice versa. The action also con- 
tains a "charging energy" term, describing phase diffusion. 
In the quasiequilibrium regime, phase coherence between the 
two states can be destroyed by either temperature or quantum 
phase fluctuations. We explicitly derive the crossover line, 
which is a function of the number of particles in the system 
A^. In the limit A^ — > oo, the internal DC Josephson effect 
is dominant and global phase coherence is established in the 
system. 

We also discuss a non-equilibrium situation by considering 
a non-adiabatic sweep of the magnetic field across the reso- 
nance. We predict that under certain circumstances, the sys- 
tem may exhibit Josephson AC oscillations. In this case, the 
numbers of molecular bosons and Cooper pairs oscillate in 
real time, slowly relaxing to a stationary state, which depends 
on the magnetic field sweeping rate. In addition, we discuss 
relevant time-scales at which such a relaxation takes place. 

Effective action. Let us consider a system of Fermi atoms 
with two hyperfine states labeled with the index cr =t> i- The 
effective Hamiltonian for the system in the vicinity of the Fes- 
hbach resonance has the following form lisll : 
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where ip and are field operators for fermions and molecu- 
lar bosons respectively, U is the fermion-fermion interaction 
strength, 6 ~ Ifi^ {B — Bq) is the detuning from the resonance, 
g is the Feshbach coupling between molecular bosons and 
fermions (which is connected with the width of the resonance 
AB d), and m and M are fermionic and molecular masses 
respectively. 

The grand partition function for the system is 



Zc =Trexp[-/3 (■-//- ;uA?)]. 



(2) 



where the trace is taken over the fermionic and bosonic 
degrees of freedom. Now let us introduce the Hubbard- 
Stratonovich field A to decouple the four-fermion term Hu 
in the action, which can be re-written in the following form 



Tt exp 



p 

-J dTHu(T) = J D^A(x)e 
J 

■ J dx[A(x)il/l- 



[-U-' fdx\A(xt] 



xTj exp 



(x)\p_Jx) + h.c. 



,(3) 



where we introduced x - (t, r) for the sake of brevity. Let us 
note that the last term in Eq. Q is very similar to the Feshbach 
coupling term in Eq. Q. To proceed further, let us evaluate 
the Gaussian inte gral over the fermionic degrees of freedom in 
the spirit ofRef. LJl. We use the following standard Nambu 
notations: 



The Hubbard-Stratonovich and molecular fields read 



A(x) 
AHx) 



and ^(x) = 



0(x) 
fix) 



The Green's function G{x) - - (t^ iAWiA' (■^)) is also a 2 x 2 
matrix. In the Nambu notations, the grand partition function 
takes the following form: 



D^A(x) D^cf 



[x) exp ! -S 



>eff 



[A,#l), 



(4) 



where the efi'ective action can be written in the compact form: 
5eff fA, ^1 = J Jr^f (r) + ^jdx \A(x)f 



-Tr[lnGo' 



-i-lnG 



(5) 



Where the trace is over the spatial and time variables and 
Nambu indices and "K^"^ is the Hamiltonian of the system of 
molecular bosons (without the Feshbach coupling). The ma- 
trix Green's function (which quite generally is a functional of 
the pairing and bosonic fields) is the solution of the following 
equation: 

■Id] 0-3- A(x) - g(f>(x) 



■I' 



2m 



G{x,x') 
^6(t-t')6(r-r'), (6) 



where is the Pauli matrix in the Nambu space. 

In what follows, we consider a homogeneous order parame- 
ter This approximation is exact in the zero-dimensional case, 
/. e. when the size of the system is smaller than the coherence 
length, otherwise one should study Ginzburg-Landau equa- 
tions taking into account the spatial dependence of the fields. 

Both the pairing field and the Bose-field are complex func- 
tions of the imaginary time. Let us write them in the form: 



A(r) = Aoe'^'<^* and 0(t) = 0oe' 



W2(T) 



and perform the following gauge transformation 







In this gauge, the Green's function reads (although the coor- 
dinate dependences of the order parameters are neglected, the 
Green's function depends on r and r'): 

xd (r, t'; r - r') = 6(t - t') 6{r - r'). (7) 

Where the new bosonic field (the order parameter describing 
a superfluid phase of molecular bosons) is given by 





Ax\-Xi) 







Now, let us consider the case when the bosonic and dy- 
namic terms in the Green's function are small compared to the 
other terms. In this limit, the effective action can be expanded 
with respect to the small terms: 

5eff [aJ] = -TrflnGo' - InGc^ + ^A^ 



+ j dT^fHr) + Tr 



ran 



-Tr 



'^2 57^ 



+ ..., 



(8) 



where Ggn is the Gor'kov-Nambu Green's function matrix: 
^ /G(T,r) F(T,r) \ 

In the momentum representation, the components of the ma- 
trix have the following well-known form I.19J : 



G(e„,p) 



'Al+€f, + s^ 



and F(e„,p) 



Al+a+el 



(10) 

where e„ = (2n + l)nT is the fermionic Matsubara frequency. 
The first two terms in action (O pin the mean-field value of 
the BCS order parameter, while the third term controls the 
magnitude of the superfluid molecular field [in principle one 
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can include higher order terms with respect to in the initial 
Hamiltonian Q, describing interactions between molecules; 
the tunneling terms in the action will remain the same in 
this case]. The most interesting contributions come from the 
fourth and fifth terms in Eq. which describe quantum 
phase dynamics. Using Eqs. and JlOl l. we obtain the 

first order correction to the mean-field action: 



6Si = Tr 



OT 



P 

= -y J dT cos \xi ir)-xi (t)], 



with 



3 gMT)^oiT) 
J(T) = -N ^^ ]n 
2 ep 



max {Ao(r), T] 



(11) 



(12) 



where is the number of particles in the Fermi system, ep 
is the Fermi energy, and e is the BCS high-energy cut-off (in 
our case e ~ ep). The leading contribution coming from the 
second order correction term [the last term in Eq. gener- 
ally has a complicated non-local form. However, if the phase 
dynamics is slow enough at the time-scale Aq', the action can 
be reduced to the familiar local form: 



6S7 



8£c J I Ot 



where 



, 3c(T) N 

/ ep 



(13) 



(14) 



and the function c(T) is given by 



c{T) 



CO 



2T 



cosh 



-2 



2T 



(15) 

with E - -yjA^^iT) + This function has the following 
asymptotic behavior at low and high temperatures: c(T «; 
Aq) = 1 and c{T » Aq) ~ 7tA^)(T)/T. Summarizing, let us 
present the main technical result of the present work: 



6S 



p 

-f 



dr 



-ycos[^,(.)-^.(.)].^(|i 



, (16) 



where the parameters J and Ec are given by Eqs. il2\ and 
J14> . The first term in action il6i can be easily recognized 
as a Josephson term, which describes coherent tunneling of 
the molecular bosons into the condensate of fermionic Cooper 
pairs and vice versa. Unlike in the conventional Josephson 
effect 1 20], the two quantum states are not separated by any 
physical energy barrier being mixed up in the real space 1 17]. 
The second term in Eq. ( II6I 1 is the effective "charging energy." 
This term describes quantum phase fluctuations. In a neutral 
superfluid, phase dynamics is possible only due to the finite 



size of the system (this is different from the case of a su- 
perconducting junction, when quantum phase fluctuations are 
possible in any system of finite capacitance). 

Quasiequilibrium regime. Let us consider first a system of 
Cooper pairs and preformed molecular bosons in the thermo- 
dynamic equilibrium, /. e. ^ - 0. In this state, the molecular 
order parameter is directly related to the BCS order parame- 
ter f21'|. The corresponding relation follows from the condi- 
tion =OandEq. Q: 



(17) 
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FIG. 1 : This diagram describes a crossover between the phase co- 
herent losephson dynamics and the phase oscillating (Pock) regime. 
The curve shows the critical number of particles as a function of tem- 
perature at which the crossover occurs. 

In the quasiequilibrium regime, the internal Josephson ef- 
fect tends to establish phase coherence between the two con- 
densates. However, even at zero temperature, this coherence 
can be destroyed by quantum phase fluctuations due to a finite 
mass term in the action ( I16> . A crossover between the semi- 
classical Josephson dynamics and quantum phase diffusion 
occurs when |22] J ~ Eq. From this condition and Eqs. J12l i. 
il4l . and illl (at zero temperature), we derive the following 
expression for the number of particles in the Fermi liquid cor- 
responding to the quantum crossover point: 



A^c(O) 



vpt/2 (6-2^) 



1/2 



,1/(vf£/) 



(18) 



where vp is the density of states at the Fermi level and we 
have used the BCS formula for the transition temperature with 
ep as a high-energy cut-off. We emphasize that the effect can 
not be called a quantum phase transition, since the transition 
is driven by finite size of the system and disappears in the 
thermodynamic limit (the system is always in the phase co- 
herent state if = 00). At finite temperatures, the Josephson 
effect is suppressed by thermal fluctuations. The correspond- 
ing crossover between the coherent state of two condensates 
and random phase rotation can be estimated from the relation 
J(T) ~ Ec(T), using Eqs. O and (EJ: 



(19) 
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where the function c(T) is given by Eq. il5\ and Ao(r) is the 
temperature dependent BCS gap 1 19]. The A^c(7^)-dependence 
was evaluated numerically and the corresponding curve is 
shown in Fig. 1. 

N on- equilibrium dynamics. Let us consider a non- 
equilibrium case, when a magnetic field is swept suddenly 
from the BEC to the BCS side of the resonance. In this case 
one has a non-equilibrium situation: the distribution function 

of quasiparticles /(k) is not given by the familiar equilib- 

-1 



rium distribution /o(k) 



The tech- 



H-exp(/J^ 

nique used previously is formally valid only if the time-scale 
at which the order parameter changes is much larger than the 
quasiparticle relaxation rate 3> t^. This is definitely true in 
the close vicinity of the transition point {A^Tc ~ A(r = 0)]: 
Ta ~ TeTc/A » Tg, where rj' is the quasiparticle relax- 
ation rate, which in a clean Fermi liquid can be estimated as 
Tj' ~ vpUT^/ep. However, away from the transition point, 
the domain of applicability of the theory presented earlier is 
determined by the following condition: A «: vpUT^/ep. At 
very low temperatures, the quasiparticle distribution function 
cannot follow the changes of the pairing field; to describe such 
a non-equilibrium regime, it is necessary to use the Keldysh 
technique 112 3ll or study the equations of motion for the super- 
fluid phase (Josephson equations) coupled to the kinetic equa- 
tions for the one particle distribution function. Unfortunately, 
the corresponding derivation is quite cumbersome (and will 
be published elsewhere). Below, we present only qualitative 
results discussing several non-equilibrium regimes possible in 
the system of interest: 

It is natural to assume that by sweeping the magnetic field 
across the resonance, we elevate the chemical potential of the 
bosonic component relative to the chemical potential of the 
BCS condensate |24]. The corresponding chemical potential 
difference 6 yields an effect similar to the effect of voltage 
in superconducting junctions. Therefore, if 5 2A the sys- 
tem should exhibit oscillations associated with an internal AC 
Josephson effect: 



Nb = J sin (6t + 0o). 



(20) 



These oscillations correspond to coherent tunneling back and 
forth between the two quantum states (condensates of molecu- 
lar bosons and Cooper pairs). Due to the AC Josephson effect, 
the number of molecular bosons would oscillate in real time 
and this behavior should in principle be observable in exper- 
iment ||3|] with the help of the radio-frequency spectroscopy 
technique (l|. 

A different physical picture emerges when we consider the 
regime of a highly excited molecular state with 6 » 2A. 
In this limit, the molecular bosons quickly dissociate into 
one-particle excitations and the Josephson effect disappears. 
In this case, one would expect to observe Rabi oscillations, 
which were described in details in Ref. 1 25] for a generic BCS 
problem. 

The internal AC Josephson dynamics described above and 
the Rabi oscillations of Ref. ]25] are dissipationless, since 



Nb(t) 




exp[-t/Trei] 



t 



FIG. 2: Real time dynamics of the number of molecular bosons in 
the Josephson regime: d « 2A. The solid curve shows AC Josephson 
oscillations in the number of particles. The oscillations decay with 
time and the system tends to an equilibrium state, which generally 
depends upon the magnetic field sweeping rate. 

equilibration is possible only in the presence of relaxation pro- 
cesses in the system of one-particle excitations. If the temper- 
ature is low enough (keeping in mind that the word "tempera- 
ture" has somewhat ambiguous meaning in a non-equilibrium 
case) the relaxation to a final equilibrium state is quite slow 
(see Fig. 2) and the corresponding relaxation rate is likely 
to be exponentially small: ttJ ~ r^' exp [-min {(5, A) /Ti], 
where Ti is the initial temperature before the magnetic field 
sweep. 

Another important question concerns the final equilibrium 
state to which the system relaxes at f — > oo. This state de- 
pends on the final temperature Tf of the system. The latter is 
related to the details of the magnetic field sweep. By rapidly 
(nonadiabatically) changin g the magnetic field, one puts the 
system in an excited state f25', '26']; in the process of relax- 
ation, the system heats up being unable to drain the extra en- 
ergy. The temperature change relates to the energy pumping 
rate 6 ~ (dB/dt) {23^2^. If the final temperature is higher 
than the superconducting transition temperature Tf > T^, the 
equilibrium phase of the system is a normal Fermi liquid state; 
in the opposite limit < Tf, the system relaxes toward a qua- 
sistationary state, which is either the Josephson coherent state 
[if the system is large enough, > NdTf); see Eq. il9\ and 
Fig. 1] or the Fock state [if < N^Tf)]. 
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